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I. tntroductjon 

One of the missions of the Aerospace Medical Research Laboratory (AMRL) 
is to define the effects of stress encountered in Air Force missions upon man, and to use 
this understanding to develop means to alleviate the harmful effects of the stresses or to 
advise aircraft designers of the limits of man's capacity to withstand these stresses. A 
second mission of the Laboratory is to develop an under s tan cfing of how a man achieves 
control of a complex, dynamic system, such as an aircraft, in order that aircraft designers 
ase this understanding to build work space environments, display system. , and control 
s rf.ems which optimize man's control capability. 

Until recently, accomplishment of both of these mission goals has been achieved 
through the classic disciplines of environmental medicine and physiology on the one hand 
and performance and engineering psychjlogy on the other hand. Both of these approaches 
have contributed a substantial amount of information and data to the aerospace vehicle 
designer and have also advanced the state of our understanding of man as a subsystem of 
a larger, complex system. Neither of these approaches, however, provides data or infor- 
mation concerning man as a subsystem of a larger system which can be directly equated 
to the t pe of information which an engineer uses In the design and development of an 
airframe. 

To design a control system for an aircraft, an engineer needs to know the per- 
formance specifications of the aircraft, the variables available to him for control of the 
airframe, the equations of motion to the aircraft and the associated stability coefficients, 
and the characteristics of the signals which will be used to provide inputs to the controller. 
For example, consider building a control system with one or m*»r e of these pieces of in- 
formation unspecified or unavailable. Consider forth he poe ,u>Uity of not having accurate 
Information aboift the effects of the flying environment u* m the integrity or stability of the 
components which must be used in building the controller or the components of the air- 
frame. 

Man, as a subsystem in a largei system, is usual’y a controller. He uses in- 
formation available to him to effect changes in the controls available to him in order to 
achieve some level of performance of the whole system. He must cope with both dynamic 
changes in tht, characteristics of the controlled system due to environmental effects and 
with changes in himself Aje to the effects of the environment and the motion of his vehicle. 
The primary hypothesis recently advanced by the AMRL is that adequate descriptions of 
man as a controller, expressed in equations ucable by the design eigineer, can be used to 
understand bow man is affected by stress and how man achieves control of complex, dynamic 
systems. Further, these equations can be used to optimize the performance capability of 
the entire system by lessening the stress effects and creating a more efficient man-machine 
division of labor. 
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In an attempt to describe man as a subsystem in manned we anon avstAmo 
several approaches to the modeling of the human operator have been tried. This naner 
*? parf ° rraod ** Systems Management Divisto^ 

Un " e f erles ana, y al8 “ 8 t«>l *»r modeling the human 

l “rf Ca *f 8 h ‘J* * Ws technif l ue c«» be utilized for determining the 
variation of the human transfer function under various levels of stress The AmwZi, * 0 

“ doe ® “X require the postulation of a model and then checking * 

Performance, ^tead, the method dete^hS rte htZT 
operator a model baaed on actual input and output data from a tracking experiment. 

II. HUMAN OPERATOR PERFORMANCE MODELING 


a) Previous Human ModgUn 


mathematical models exist to describe human operator oer- 
formancei A hough a study of these models leads one to conclude that thev indt at 
problems .relatively narrow viewpoints, prevtoWiwcte?^ 

knowledge and greater understanding of this complex emblem * %!£?**** 


b) AMRL Mode lh 




of multiple discrete and corttniuToStrol patte ‘ '£ h the coroptex mlxtare 
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c) Effects of Stress on the Human Controller 


30 a 


Stress fan ill-defined word but a usable concept) can affect man physically, 
such a* by acceleration forces; physiologically, such as by altered central nervous 
system function in response to drugs; and psychologically, by anxiety and fe^r. The 
modeling of human transfer functions which can account for stress Is of great interest* 
physical stresses can have a direct effect upon manual control capability. Physiologic 
and psychologic stress, however, can indirectly affect man's control capability. The 
manifestation of these indirect effects in changes in the human operator can be extremely 
subtle. There are tests, which can be given to a stressed man, which ate more sensitive 
in exposing altered human performance due to the stress than can be currently developed 
by the manual control community. Rather than being an affront to the methodologies 
used by the manual control community, this fact Is an accolade to man's capability to 
compensate for the indirect effects of the stresses. Total system performance may not 
measurably change until the physiologic and psychologic structure of the roan cannot sus- 
tain intelligent responses on his part. 

That man, as a controller, is changing wider mild physiologic stress Is assumed 
to be true; and the fault with exist teg methodologies for human operator identification is 
assumed to be their lack of sensitivity in elucidating 'hese changes. The study which is 
reported upon here was designed to expose the veracity of the above statement. 

d) The Experiment Analyzed 


A MR L has the facilities to expose man to a number of differing kinds of stress 
and combinations of stress. Acceleration stress, as one of Its effects upon man, causes 
both local tissue and system c hypoxia. In order to separate the physical effects of 
acceleration stress from the hypoxic effects, a short study was run using three subjects 
and breathing mixtures representative of aea level and 15,000 feet altitude 
(PA 02 =* 60* saturation). The subjects were instructed to perform a single axis com- 
pensatory tracking task consisting of a random forcing function and a third order repre- 
sentation r* the pitch axis dynamics of a high performance aircraft. Randomization of 
the sequence of altitude exposures and multiple replications of the runs were accomplished 
to maximize the accuracy of the measurement of the error variance. The study concluded 
with the observation that the variance and mean of tracking error for the two minute task 
was not affected by the changes In altitude. The data acquired during the two minute runs 
consisted of the forcing function, the tracking error, and the control stick output sampled 
at 20 msec intervals. This data was subsequently analyzed using time series analysis. 

Figrre 1 illustrates the experimental compensatory tracking configuration utilized. 

The forcing function, n(t), represented the resultant outptd of a noise generator. The 
transfer function of the plant, H, was set at 1. I/s, and l/s 2 . The human operator was attempt- 
ing to track the input of a compensatory tracking display in a simulated F-4 cockpit . 

For purposes of this report, the three subjects utilized for this ftudy will be 
referred to as subjects A, B, and C. In order to gain flexibility In analysing the experi- 
ments, several runs were repeated the same day and on other days in order to dete.' nb 
variations due time, learning, adaptability, etc. 
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Figure X. Experiment Block Diagram 


m* DATA REDUCTION USING TIME SERIES ANALYSIS 
a) Time Series Representation 

Time series representation of discrete linear random processes consists of 
linear difference operations which relate the time series data x t to a white nefee pro- 
cess, z fc . A purely white noise process is represe nt ed by 

V*t a» 


An autoregressive (AR) model of order p is represented by 

X t =# l X t-l + *2 *t- 2 * *1 i + **• * *p X t-p + e t (2) 

where * t for i = 1,..., p are the model coefficients. By defining the backward shift operator 


Bz 18 z 
t t-1 


then the general AR equation can be given by 

(1- B - # 2 B 2 - ••• - # t B l - ... - B**) Xj = z t 


( 3 ) 


(4) 
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A moving average (MA) process of order q is represented by 


x t 3 Z t " # 1 Vl “ *2 z t - : 


9 i n - l " * *q t -q 


In terms of the backward shift operator B, this equation is given by 
•jB 1 - ••• - # q B q > z t 

A genera! mixed AR/MA process is given by 

V Vt-l * *2^-2* — * *p*t-p ,Z t 

* *1 Z t - 1 * *2 Z t - 2 * •** ‘ *q *t - q 
In terms nf the B operator, this equation fie given by 


(5) 

( 6 / 

(7) 


(!-♦, B- * 2 B 2 b‘ ♦ p B P ) ^ = (1- • 1 B-* 2 B 2 


B h )z. 


Differencing of the data is necessary to induce statlonartty, if the data Is 
comprised of nonstatfonary processes such as a ramp or random walk in addition to 
a stationary' process. The representation of one difference ot data Is given by 


(1 - B) x t 

The represents lio-i of d differences ot data is given by 

< 1 - B> d x t 


(9) 

( 10 ) 


b) Model Building Approach 

Data reduction using time series analysis prefaces a mathematical model by 
reducing a wave form to white noise while identifying the correlated portion of the time 
series. The mathematical model is obtained by s ree stage iterative procedure based 
on identif icat Ion, estimation, and diagnostic checking. The identification process is 
concerned with the generation of the series in order to determine a class of models that 
should be investigated. The estimation phase uses the data to make inferences about 
parameters conditioned on the sufficiency of the model chosen. Diagnostic checking 
anaLyr *s the fitted mode 1 with the data in order to determine any model Inadequacies 
and obtaining model Improvements. 


c) Identification 


In the identification stage of model building, the observations x. are differenced 
as man;, times as required to induce stationer iy. The resulting AR/MA process is then 
identified by analyzing the autocovartance function (ACVF), the autocorrelation function 
CACF>, and the partial correlation function (PCF^M These functions are also useful for 
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obtaining approximate estimates of the parameters which are useful in the estimation stage 
for providing initial values needed by the iterative estimation procedures. As sh jwn in 
Figure 2, the pure AR and MA processes are duals of one another, and opposite behavior 
of ACF and **CF indicates a pure AR or MA process. 

d) Estimation 

The estimation technique for determining the# and 0 parameters used Marquardt's 
Algorithm 131 . This nonlinear estimation approach performs an optimum interpolation between 
two traditional nonlinear estimation techniques - the Gauss (Taylor Series ) method the 
method of steepest descent. The algorithm combines the advantages of the ability of the 
steepest descent to converge from an initial guess which may be far from the final value, 
and die ability of the Gauss method to close in cm the converged values quickly once the 
vicinity of the final value has been reached, in addition, Marquardt's algorithm overcomes 
the disadvantages of the slow final convergence of the steepest descent and the possible 
divergence of the successive iteration with the Gauss method. 

e) Diagnostic Checking 

The diagnostic checking step Is concerned with examining the white noise residuals, 
zj, from the fitted models in order to determine any lack of randonmess. A nal ysis of i* 
can incficate whether the model is adeqr or inadequate, and yield information on how to 
better describe the series. Modeled mourns would then be refitted and residbfected to 
diagnostic checking. 


AS PROCESS UA PROCESS 



Figure 2. Comparison of ACF and PCF for Pure AR and MA Processes 
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Transfer Function Development 


Time series analysis techniques can be extended to obtain d iso ret e linear trans- 
fer functions of systems having an Input time series x and an output time seriee y t , For 
th^. system Illustrated in Figure 3 t it 'n assumed that the system can be adequately 
represented by the impulse response parameters, v , v, , v of v Q , v*, and v., There- 
fore, y^ is related tox^by oiase o 

y t = <v o * V 1 B + v 2 ® 2 * v 3 ® 3 * v 4 ® 4 + v 5 B ^ x t fll) 

In general, for a dynamic system containing input observations x. and output observations 

v . • 


y t = V<B), H (12) 

where v (B) represents the impulse response function of the system. 

The Input/output relationship of the type of system represented by Figure 3 can 
also be represented by the following general linear difference equation! 

"-fi® f H® R >y. *<“o* W l® (13) 

f (B)V“(B»V b (14) 

Solving for y t , the following expression Is obtained: 

y t *>(B) (B) 01 (B) B b x t (16) 



Figure 3, The Input/Output Model 
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a«> 

(17) 


tiling the dnflnltlon of Equation 12 tor the ££5* *"* 

and substituting this tote the difference Equation 13, the follow tog identity le outtlnedi 

(1- f, B B®) (tr, * VjB+VjB®*...) - "B® 8 ’ <U) 


The following set of equations Is obtained by eqputog Ifco power* of Bt 


». J<b 


at) 

f l V j-l * f 2 V i-2 + “* 

f R V J - R *% • » * b 

(% 0 ) 

f l v j-l * »**»-**•” 

* !R v J-R-* , J-b' , ' b * 1 ' ”• 

, b ♦ 0 (21) 

f l v J- 1 * f 2 V J-2*'" 

• T R Vr .)>brS 

(22) 


The retatlonehlp of the v peremetera from the impuleo r oaponoo £ £££ £1 X* * 
parameter! of the dUferino* equatten rapraaenutlon on Mint rated ieE*i Mona 12 
through 22. 

The traMfhr taction model Identification fiMNi *• ilmptifted tfthe 

input to the eyetom le white noiee eiooe the oreea oorwi mta tao^ bmw nnn ». 
nolee input nod me eyetem output rapreeeota the lnptd— 

V the input le com* ether etoobaeMo progeM, al mplHleWeo l e B»n«*w by jwymang 
an operation denoted a* "iprewhttenlng".'*' TMe eppreaeb wea need to *ln atedy. 

The eatimettoo teoholque need for obtatotng the f •odotpernmetera tadWwwmo 

Equettcol2laaonlo*auatetba20v«f*rt|m«dlmeoitaeln»odittof oaee. AfUer ^ptylnc 

Marquardt'a Algorithm, the reaultlug medal le a* the term 

y,» f " l <S) «<*> B b Xj ♦ R t (22) 


whare R. repreaenti the raatdual Ugoal not aorreleled with the Input. Sobrtng eqmtlor 
22 for r|, the following axpreealen fa obtained! 

R t • [ y t - f " l (B) w (Bt B b x, ] <*♦) 


S» 




lum 'zgzissr ,h * fand,upar *' ne<, ’ r ‘ * *» .»» 

s(r.w) j r, 2 

I *| 1 <2fi> 

l^HiagnoBtic checking procedure is analogous to that given tor the one dimensional 

K‘ Transfer Function Development In a Feedback Loop 

ThU study was concerned with modeling the human operator in a compensatory 
tracking system ■ as shown in Figure 1. The human transfer function 0 cannot be ob- 
tained directly by only considering th, time series et and u t , since the residual Rt muat 
be independently distributed of For example. Ha orosa correlation function muat be 
zero for all lags: 


f. [fR, -R) (c |hT -c) ] o 


If this is violated, then the parameter estimates obtained are not consistent. Unfortunately. 

3 9yate "\?» **' * bum< ”> tr»n*t*r function, is Joosted within a feedback loopf 
a ! 8Umpt ^ n ®t »nd R t does not hold. Ml In order to overcome the 

problem, the signals et and il must be analysed with respect to the Input signal R* Since 

ritdwt e Kf e r 1y f ner !!f drandom It Is^’orrehTted^th gone- 

rattd within the system and present at the output or error signal. 

relating *>"<>"•»« «n»tlon. c.n b. define 


"» ' ♦ 6(8) H(B) "t * *1, t 


*» 1 * C(B> H(B) "t *2, t 


G ( 8> Human transfer function 
H (B> Plant transfer function 
I t* z 2 t No,s * components on the output 

signals uncorrelated with the Input signal n. 


B> redefining the control ratios 

T » ,r * r^WTHBi 


T * <B > ■ T^nw H(B) 


(30) 


«>. human tranafor function 0(B> o.n b. obtain* u.lng the following aquation. 



*" ° W * in * d fro,n ‘ ra,uf,r 

IV. HUMAN OPERATOR PERFORMANCE MODELS 
r..-.,., th. MM.* JLt mJL can biS 

* danotn'lnator. 0 ** * Mm * Mmr * °' 1 "•«* " ,th * ”">"<1 ord.r 

• Two of tho three models hod a aero in the numerator. 

# *%SgSl£'?* m ^ —— — W in ovar- 

& ST * *••• -W. wm that 

clualon raaohad on thl. ah* which 5 

ralatad JSSKSS&STS^ * * ’**'**" »'«*• - 
ate the seasonality character let ioe from the dua manner* In order to ellmin- 

the data by the pertoT^ T iify 10 

and emphaeiaee the actual reepeo^T he refers lh# ••••°n*Uty pattern 

ticulariy important role to the analysts of the B h * Vi plBytd • P« p - 

that the simplifying operattoTp^ -1 * f.ft ^ ^ 

au^gS^g ^ow saasaas 

r.££. tTm 

raaponaa In tracking I— »bp»lmpa.ln t on thalr 

investigated further. ^ m * tr,,i itwc * and its origin needs to be 

punt'. wi " d * — ■ «*-"«•«*« •- wm* «. 

H •> -L 

' (3*1 


hs- ■*”' 
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in order to dewrmlno Iho change in the human ijrawta* ifiEdlia 

In order to compare j* Jormanc*' the oxperimont the! aub)eet C performod to , a!! 0 ’ 
00014, Run No. 00001, was acalytod and ooinpered with Ma JL 

Record No, 0J004, Run No, 00* rt , uatof the time 

modela obtained are compered to iguro 6, Anelyola of tbeae two model* hmtoeM that the 
moat atgnifl >*nt difference whh H - 1/e e^nepojtntlon 

numerator* loading one to conclude that the toteg rotor, H • 1/a, prow Idea memory to the 

•***•«• tse s 

the two numerator lead factor* of (1-0,B62B) and (1 ♦8.76W, However , etoce the ooeffi 

oiecN of the B ten*; 0,7ft, l* very much greater then the coefficient of 1 In the term 

a + 8?7KB) , It to Interpreted to^ean thitalthough there le a .mall toatontoooeua reepenee, 
the predominant reaponae le delayed by 0,2 tee. There tore, the numerator doee appro- 
ximate a transportation log with the following reeutti 


- 3, 02 ( 1 - 0 

fl-MWW 


m 


(33) 


The Interpretation of the tocreeee in operator gain from ^ 

lion with H ■ 1, to 3.08 In tho appronlmatod MWjrmiltr tootloawIthH - U* 

to the deoreoee In open-loop gain oaueed by the H ■ 1/a, Thoe, when the operator la coo- 

^m^TplMn H Tm/.Tu to nacaaaary ttot to tonu. hto ,.to to t^ no.to tor tto 

^i puaHMi Mutod by tto plant In ordor to aoblava • daalrabla bandwidth. Thto proM«n 

waa not preaant with H - 1. Tharatora, Hto Mat tto otong* n tto p Into f raw »M 

to 1/a aupptlaa tha oparator with memory from wbtob to la tottor aMa to anttotpota i Mura 

notion, and alto daoroaaoa tto opm-loop tain oauatog tto oparator to Inoraaia hto gain in 
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Figure 6. Seasonality Periode Oenereted by the Hu men Operator 
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Figure 6, Model* of Subject C with Plante uf H-l and 1/a 


-53* 


-334- 




analyzing expe rtmenu° a i m*| a r'uM toi o cone UteJ£d'i *“*«*«> waa atudlad by 

subjects A, B, gnd c related ZT. ^ P1 « ur “ 4 *nd S, where the three 
applied to examine 30-secood ’Vtndewe^JaM^lim «Da?ad**| , n' , Th * pro * rRm w » < alao 
the duration of the three ordinal experlmUte wlya 2 T£**£J ; 3 l. ,n ** rV * U> ov#r 

; • “ubj J,° A^; r nt! ,h *. ,r *'“ , * r , ' wct '°"® 

C repeating Record No. 00004. Run No. 00001 Imn '«<il«tely, and for eubject 

7 compares (he two models obtained for »e?uV» I ,# !SL 4,y ! *" Up w * re Figure 

obtained (or subject c. ODWItoed for •*»«>* A, Figure 8 compares the two models 
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Figure 7. Repeatability of Subject A with H ■ 1 
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by haying'a^j^V^mme^iai^r^^^'^V'^®’^ "**« A he. changed conaidarably 
transportation lag does not *«(■! » l. ^perlment, ft it signlfioant to note the ■ * 
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the human operator fa applying much more smoothing during the earlier part of the ex- 
periment than at the end. Therefore, hie reeponae time fa faater towarde the later part 
of the experiment indicating that he haa learned and can anticipate change*, with the re- 
sult that he can reapond faater. 

Time series analysis haa the capability of modeling human operator transfer 
functions under task and/or environmental stress conditions. At an illustration of its 
application In this manner, consider a human operator In a compensatory tracking situ- 
ation, It is desired to detormina the variations of his model’* parameters under hypoxia 
stress and he ia required to perform the trucking task under simulated oxygen leva's 
corresponding to aea level, 12, 000 feet and 22,000 feet. For each of thess three exper- 
iments, the program can determine the human transfer function directly from the dat* 
vary eff^lcr/Jv, 

It might also be desired to determine the parameter variation of the human trans- 
fer function under combinations of stress. For example, the program could be extended 
and applied for determining the human operator model while ha la tracking a target on z 
compensatory display with various combinations of stress Including hypoxia, roll, pitch, 
vibration, heat, etc. The program could be further applied for determining the effect of 
these stresses on such physiological parameters as blood pressure, heart rata, etc. 

Due to budget limitations the program was not extended to model stress data. 

For the srtme reason, the program was not applied to model pilot compensatory tracking 
data for areas where the element H was changed to l/s 2 . 

CONCLUSIONS 

This pioneering effort for applying time series analysis to model the human 
operator represents the first time that time series theory 4 has been applied to this problem 
and Its results have been most Illuminating, The foremost conclusion found on this study 
Is that the time aeries method is a very efficient, effective, and powerful method for 
modeling any dynamical process having an iignit and output which contains noisy observa- 
tions, Application of the time aeriea technique to modeling the human operator in s 
compensatory tracking situation Indicated that he can be adequately represented In the B 
domain by the B- transfer function having the following form. 



m) 


plus a rosldual containing s completely modeled mixed AR/MA structure and white noise. 

It wts also concluded that when the plant's transfer function was changed from 
1 to 1/e, the integrator, H 1/a, provided memory co the system, and the operator was 
better able to anticipate and predict future actions. The resultant modal, therefore, 
did not contain a transportation lag. However, a good approximation to this transfer 
function contained a tranportatlon lag and it was found that the operator also had to In- 
crease his gain significantly. The Interpretation of his Increase In gain was caused by 
hie compensation for the decrease In the open-loop gain due to the integrator H * 1/s, 

The operator basically increased hla gain In order to achieve a desirable bandwidth. 

It was also found that when an operator Immediately repeated an experiment, he 
attempted to eliminate the transportation lag and one of the smoothing time constants. 
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This is a very significant result and is anticipated baaed on his ability to learn. However, 
due to the magnitude of the coefficient, the model atilt iad an effective transportation lag 
from a practical viewpoint, ft was also found that having an operator repeating an experiment 
several days later Indicated very tittle learning ability. 


Analysis of the variability of the human transfer function over the length of the 
experiment by means of analyzing a small number of points (window) indicated that the 
human operator Is applying much less smoothing during the later part of the experi- 
ment. This was interpreted to mean that he has teamed and can anticipate changes, with 
the result that he can respond faster. 

A very significant result of this effort was measurements of a seasonality being 
generated by the operator. These seasonalities are interpreted as s rhymlc effect that 
the operator is generating in order to trade random noise. 
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